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We analyze exciton spin dynamics in GaSe under nonresonant circularly polarized optical pumping 
with an exciton spin-flip rate-equation model. The model reproduces polarized time-dependent 
photoluminescence measurements in which the initial circular polarization approaches unity even 
when pumping with 0.15 eV excess energy. At T = 10 K, the exciton spin relaxation exhibits a 
biexponential decay with a sub-20 ps and a >500 ps time constants, which are also reproduced by 
the rate-equation model assuming distinct spin-relaxation rates for hot (nonequilibrium) and cold 
band-edge excitons. 


I. INTRODUCTION 

High spin polarization and long spin-relaxation times 
in solid-state systems are desirable for many spin- 
tronic applications. Toward this goal, nonequilib¬ 
rium spin dynamics have been investigated in various 
semiconductors^^®, with gallium arsenide (GaAs) the 
most studied example. However, the degenerate heavy- 
and light-hole valence bands and sub-ps hole-spin re¬ 
laxation in bulk GaAs have limited optically pumped 
electron-spin polarization in GaAs to 1/2 and the degree 
of circular polarization of the resulting photolumines¬ 
cence (PL) to 1/4^’^. The electron-spin relaxation times 
can be enhanced by doping' or quantum confinement®. 
For example, near unity electron spin polarization can be 
obtained in heterostructures where heavy- and light-hole 
energy degeneracy is lifted by quantum confinement or 
strain®^^^. 

In a recent report^^, we demonstrated that generation 
and preservation of a high degree of optical spin polariza¬ 
tion, p, is possible in nanoscale slabs of GaSe even when 
optically pumping 0.1 to 0.2 eV above the gap. This is 
a consequence of the unique bandstructure (Fig. 1) of 
the group-HI monochalcogenides^®’^^, which consists of 
orbitally nondegenerate uppermost valence band (UVB) 
and lowermost conduction band (LGB) with minimal 
spin-orbit-induced spin splittings^®. Most importantly, 
in GaSe, the large separation of the LGB and UVB from 
adjacent bands results in the suppression of Elliot-Yafet 
spin relaxation mechanism. In this study, we show that 
the exciton spin dynamics in GaSe at low temperature 
can be quantitatively reproduced by a simple model in¬ 
cluding momentum relaxation of excitons and spin relax¬ 
ation of excitons and carriers. 


II. CRYSTAL AND BAND STRUCTURES 

The HI-VI semiconducting compounds GaS, GaSe, 
and GaTe (MX) all form layered crystals. The bonding 
between the layers is weak, resulting in the easy cleavage 
of these crystals. In the case of GaS and GaSe, each layer 
consists of four planes of atoms in the sequence X-M-M- 
X and belong to the space group Different stacking 


orders of hexagonal layers result in formation of four com¬ 
monly known polytypes: (3 (Dg^), rhombohedral 

7(C'3„), and J(C'I^). Symmetry-dependent optical transi¬ 
tions and spin dynamics near the band edge can be sensi¬ 
tive to polytypes^®. Here, we study e-GaSe which has an 
ABA (Bernal) stacking order and belongs to space group 
— P6m2, which is noncentrosymmetric, i.e., lacking 
a spatial inversion center. 

Bulk GaSe is generally regarded as an indirect-band- 
gap semiconductor. An indirect transition from the T- 
point to the M-point ^10-20 meV below the direct gap 
at the F-point is nearly resonant with the direct exci¬ 
ton transitions at the F-point (exciton binding energy 
20-30 meV)^®“^®. Fig. 1 shows a schematic crystal struc¬ 
ture and bandstructure of bulk e-GaSe. The LGB and 
UVB near the T-point are derived primarily from Ga s- 
like and Se p^-like orbitals and have respectively r4 and 
Fi symmetry. The valence bands are split by crystal- 
field anisotropy and spin-orbit interaction, leading to two 
bands with Se Px,Py symmetry about 1.2 and 1.6 eV be¬ 
low the UVB maximum. 

The strongly anisotropic crystal structure leads to cor¬ 
respondingly anisotropic optical properties. The polar¬ 
ization vectors E || c and E T c belong to the F4 and 
Fg representations, respectively. Gonsidering the tran¬ 
sition between UVB (Fi) and LGB (F4) without spin, 
the direct product Fg x F4 belongs to the representation 
r4 of D^h: therefore, the direct transition Fi —r4 is 
orbitally allowed only for E || c. Taking into account 
spin, the direct transitions then occur between valence 
and conduction bands with the following symmetries in 
the double group D^h- 

r4 —Fg (s-like conduction band) 

Fi —)• r7 {Bjp^Aike valence band) 

Fg—)'r 7-|-Fg {A,px^y-\\ke valence band) 

Fg —)■ Fg -I- Fg {C,Px,yA\ke valence band) (1) 

The direct product r7xrg = r3-|-r4-|-Fg then contains 
the representations for both E T c and E || c so that the 
optical transition between the UVB and LGB becomes 
weakly dipole-allowed for E T c. 

A complete accounting of the optical response requires 
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FIG. 1. Band structure and selection rnles. (a) Crystal 
structure of e-GaSe showing ABA stacking of the individual 
layers, space group Dlf^/P6m2 (#187) and point group Dsh- 
An individual layer consists of four planes of Se—Ga—Ga—Se, 
with the Ga—Ga bond normal to the layer plane arranged on 
a hexagonal lattice and the Se anions located in the eclipsed 
conformation when viewed along the c-axis. The inner solid 
boxes represent a unit cell, (b) Sketches of the band structure 
of e-GaSe at the F-point and the representations to which 
the states at the F-point belong with (without) spin-orbit 
interaction, (inset) Direct-gap excitons and selection rules. 


an excitonic (two-particle) picture. To a first approxi¬ 
mation, one can regard band-edge excitons derived from 
carriers only in the conduction band and the Fi UVB. 
This structure is then perturbed by spin-orbit mixing 
of the Ts valence band into the UVB and then further 
perturbed by electron-hole exchange. The direct transi¬ 
tion from the electronic ground state to an electron (Tg) 
and hole (Ty) in an s-like (Ti) direct-gap exciton state 
yields an exciton belonging to the following (degenerate) 
representations^^’^"^’^®’^°’^^: 

= Tt X Tg X Ti = r4 + Tg + Te, (2) 

where r 4 is a pure spin singlet and Tg and Tg are pure 
spin triplet states. 

Including spin-orbit mixing of the Tg lower valence 
band into the UVB, the s-type band-edge excitonic wave- 
functions become^® 

r4 = |5)+a4|T), 

Ta = \T ), 

Tg = |r)+a6|^), (3) 

where \S) and \T) represent singlet and triplet states, 
respectively. Using Ago measured by Sasaki et 


al.®^’®^, the coefficients 04 « ag « {Aso/Eba) ~ 
(0.44eV/1.27eV) « 0.35. The Tg state is mixed into 
the UVB at about the 10% level, consistent with the ra¬ 
tio of the oscillator strength (absorbance) of E T c and 
E II c. The upper level r 4 corresponds to a total exci¬ 
ton spin S' = 0 and experiences a splitting Ai « 2 meV 
because of electron-hole exchange^®. The states Tg and 
Tg correspond to a total exciton spin S = 1, and Sz= 0, 
±1 and are nearly degenerate (energy splitting A « 0). 
These states are thus labeled by the indices 0 and ±1. 
The r 4 state can be excited by light with E || c. For 
optical excitation with wave vector k || c, Tg (S^ = ±1) 
states can be excited by circularly polarized light with 
E T c, whereas the Tg state is optically inactive^^h4^ 
Earlier studies suggested that a high degree of optical 
orientation in GaSe could be achieved under nearly reso¬ 
nant excitation at low temperatures. Using steady-state 
measurements, Gamarts et al.^^h4 demonstrated optical 
orientation and alignment of excitons in GaSe by showing 
luminescence with circular polarization above 90% under 
steady-state circularly polarized optical excitation in res¬ 
onance with direct excitons at cryogenic temperatures. 


III. RESULTS 

Sample preparation and experimental measurement 
methods were described previously^®. Thin films of GaSe 
crystals were mechanically exfoliated from a Bridgman- 
grown crystal®® and deposited onto a silicon substrate 
with a 90 nm Si 02 layer, with the thickness measured 
using atomic force microscopy. Samples were mounted 
in vacuum on a copper cold finger attached to an opti¬ 
cal liquid helium flow cryostat for all experiments. GaSe 
nanoslabs were optically excited by 2 ps laser pulses from 
a synchronously pumped optical parametric oscillator 
(Ap ^ 560-595 nm, Ep ~ 2.21-2.08 eV) or by second- 
harmonic pulses from a Ti:sapphire oscillator (Ap ~ 410 
nm, Ep ^ 3.0 eV). The laser beam was focused through 
a microscope objective (numerical aperture N.A. = 0.28) 
to an area of about 80 /rm® on the sample. The wave 
vector of the pump is along the crystal c-axis (the sur¬ 
face normal), i.e., the electric field vector E is orthogonal 
to the c-axis (E T c). The polarization and the flux (P) 
of the pump laser were controlled by liquid-crystal-based 
devices without mechanical moving parts. The samples 
were excited with the pump flux P from 0.1 Pq to Pq, 
where Pq = 2.6 x 10^"^ cm“® photons per pulse. We esti¬ 
mate the photoexcitation density to be from « 2 x 10 ^® 
cm“® to 3.4 X 10^^ cm“® (2.7 x 10“^® cm“® per layer) con¬ 
sidering the absorption coefficient at 2.1 eV (« 10® cm“^ 
for E T c) and Fresnel loss from reflection. The pho- 
toexcited carrier density is below the Mott transition®^ 
of direct excitons occurring near electron-hole (e-/i) pair 
densities of 4 x 10®^ cm“®. 

Time- and polarization-resolved PL measurements al¬ 
low us to determine separately the recombination time, 
the initial spin orientation, and the spin-relaxation time. 
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Polarized PL measurements were performed under exci¬ 
tation with excess energy of about 0.1 to 0.2 eV above 
the exciton emission peak. The band-edge exciton PL 
emission at room temperature is near 620 nm (2.0 eV), 
independent of thickness. The exciton peak gradually 
red shifts from 590 nm (2.1 eV) in bulk (thickness> 1000 
nm) to 620 nm (2.0 eV) to 90-nm nanoslabs at T = 10 K. 
This shift is attributed to increasing contribution to PL 
from localized excitons, which is the subject of ongoing 
studies and beyond the scope of this paper. In addition, 
we observe that the quantum yield of luminescence is 
greatly suppressed in sub-50-nm thick samples (Fig. 2). 
In this paper, we focus on a 540 nm thick GaSe sample 
(S540nm). 



FIG. 2. PL quantum efficiency. Quantum efficiency of 
luminescence as a function of thickness of GaSe nanoslabs at 
T = 300 and 10 K. Here, quantum efficiency is defined as 
the ratio between the luminescence emission flux and opti¬ 
cal absorption flux per layer. Optical absorption flux is de¬ 
termined using experimentally measured reflectance for each 
sample and reported absorption coefficient a = 1.1 x 10^ 
Optical collection and detection efficiency is mea¬ 
sured by passing a 633-nm laser beam with known power 
through the optical set-up and spectrometer. The emission 
flux is then calculated by including Fresnel reflection loss at 
the surface assuming angularly isotropic emission. 

Fig. 3 shows time-integrated PL spectra under co- and 
cross-circularly polarized excitation and detection at T 
= 10 K. The PL is peaked at about 2.07 eV with a low- 
energy tail that is attributed to localized excitons. As 
in samples described previously^^, the degree of steady- 
state spin polarization {pc) is highest at high energies, as 
expected given the reduced time spent at such energies. 
There is also a pronounced increase in pc as the excitation 
photon energy decreases towards the band gap. This is as 
expected from the greater spin-orbit field at larger k in 
the D’yakonov-Pereh (DP) mechanism of spin relaxation 
and the increased rate of scattering at larger k in the 
Elliott-Yafet (EY) mechanism. 

For more direct insights into the nature of the spin 
relaxation, we performed time-resolved measurements of 



Energy (eV) 


FIG. 3. Polarized PL spectra under optical excita¬ 
tion at 2.156 eV and 2.214 eV. Time-integrated PL spec¬ 
tra (A)(co-circular, dashed lines) and 1'^^“ (A)(cross¬ 

circular, dotted lines)] and degree of circular polarization 
Pc{E) (solid lines) of S540nm GaSe samples under ex¬ 
citation at pump flux P = 0.5 Po, where Po = 2.6 x 10^"^ 
cm~^ per pulse. Blue (Red) lines are for Epump = 2.214 eV 
(2.156 eV). 


the spectrally integrated PL (Fig. 4). The PL at cryo¬ 
genic temperatures (10 K) reveals a fast (<10 ps) rise 
time followed by a biexponential decay with time con¬ 
stants Tg « 20-50 ps and Tq « 150-200 ps in all sam¬ 
ples. The time-dependent circular polarization pdt) is 
also observed to be biexponential with decay time con¬ 
stants r' « 30-40 ps and r" > 150-200 ps. The ini¬ 
tial, resolution-limited rise in the PL signal represents 
the rate of carrier scattering to small values of the mo¬ 
mentum (i.e., to near the band edge at the P-point). The 
biexponential nature of Pc{t) can be understood as a con¬ 
sequence of fast spin relaxation during thermalization to 
the band edge followed by slower spin relaxation once the 
carriers are at the band edge. The initial decay of spin po¬ 
larization is greater at 2.214 eV than at 2.156 eV, again 
consistent with the expectation that greater spin-orbit 
field and faster momentum scattering at higher energies 
should result in faster spin relaxation. 

Previously, we showed that carrier spin relaxation 
leads to rotation of the emission dipoles from in-plane 
to out-of-plane, i.e., from perpendicular to the c-axis to 
along the c-axis^®. This is manifested by the polarized 
remote-edge luminescence (REL) 100 pm away from the 
pumping spot at the cleaved edges of the sample as a 
result of index-guided (waveguided) emission from out- 
of-plane dipoles^®. Time-resolved measurements of the 
spectrally integrated REL at T = 10 K reveal an initial 
rise of the REL with a ~30 ps timescale that matches 
the initial fast decay of the PL from the excitation spot. 


IV. SPIN RELAXATION IN GaSe 

Spin-relaxation mechanisms in semiconductors depend 
on details of the energy versus momentum dispersion 
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FIG. 4. Polarized PL dynamics under optical exci¬ 
tation at 2.156 eV and 2.214 eV. (a) Time-dependent 
PL intensity and (dotted red, 

blue, and black, respectively)], degree of circular polariza¬ 
tion pc(t) = ^+/t|t)7r+/-(!) magenta), and R{t) - 1 = 

“ 1 (solid cyan) under excitation Epump = 
2.214 eV at P = 0.5 Pq. (b) Same as (a), but for excita¬ 
tion Epump = 2.156 eV. The ±/± [X/Y) labels represent the 
helicity (polarization) of the pump laser and luminescence, 
respectively. 


and spin-splitting because of spin-orbit interaction^®’^*^. 
In III-V semiconductors such as GaAs, the group-V 
Py, and Pi;-derived valence bands are closely spaced be¬ 
cause of the small crystal field. Consequently, the angu¬ 
lar momentum and quasimomentum of holes are strongly 
coupled, and the spin orientation of holes is lost in a 
period comparable with the momentum relaxation time 
The slower relaxation of electron spin in 
such semiconductors is a result of the reduced spin- 
orbit mixing of the group-III s-orbital-derived conduc¬ 
tion band with distant bands. Electron-spin relaxation 
is usually analyzed in terms of three mechanisms^’®’^’®': 
D’yakonov-Perel’ (DP)^’®®’®®, which is associated with 
spin-orbit-induced spin splitting of a band in noncen- 
trosymmetric systems®®; Elliot-Yafet (EY)"^®’^®, which 
is associated with spin-orbit-induced band-mixing; and 
Bir-Aronov-Pikus (BAP)^®^^®, which is associated with 
electron-hole exchange interactions. For photoexcited 


carrier densities above 10®® cm“® as studied in undoped 
GaSe here, we neglect spin relaxation due to the BAP 
mechanism. 

The EY and DP mechanisms are associated with the 
spin-orbit interaction (SOI) and the spin-orbit-induced 
spin splitting, As(k) = |if(k,t) — i?(k, 1)]®’®'’. Using a 
four-state (two bands with spin) model Hamiltonian in 
the absence of an external magnetic field, one can relate 
the spin-relaxation rate of electrons (holes) with quasi¬ 
wave vector k away from the conduction (valence) band 
edge with the following equation®'’: 


|A.(k)p rp|T(k)|® 

Pp ^r® + A2(k)’ 


where Pp = fi/rp is the scattering rate of electron/hole, 
with Tp being the corresponding momentum scattering 
(or correlation) time, As(k) being the spin-orbit-induced 
spin splitting, and L(k) being the SOI between the adja¬ 
cent bands with energy separation Ag. 

In GaSe, the initial carrier cooling to the band edge oc¬ 
curs in the sub-ps to sub-10 ps range (Fig. 4)^®. This fast 
energy and momentum relaxation (Tp ~ 1 ps) is demon¬ 
strated by the sub-10 ps PL rise time that is nearly inde¬ 
pendent of PL emission energy. The Pz-like UVB is well 
isolated from the LVB and the adjacent p^^y-lWe valence 
bands (i.e., Ag ~l-2 eV), and as a result L/Ag « 0.02- 
0.04'®®’’^®. The hole-spin relaxation caused by the EY 
mechanism Tf®^ « (^)^rp (Pp ^ Ag(k)) is expected 
to be much smaller than the momentum relaxation rate 
Pp, resulting in a spin-relaxation time of ^1000 ps or 
more for Tp ^1 ps. 

In e-GaSe, hole spin splitting Ag(fc' = 0.15) « 5 
meV®®, where k' = |k|/rAr. As a result, for holes with 
finite momentum k’, |A^(A:')p/rp rp|L(A:')p/A®(fc'), 

i.e., spin relaxation of hot holes in GaSe is dominated by 
the DP mechanism, in which spin relaxation occurs from 
the precession of spins in an effective magnetic field as¬ 
sociated with Ag (/c')®’®’®"^. When the DP mechanism is 
dominant, the smaller the spin splitting, the longer the 
spin relaxation time « h/T^ for the same momentum 
relaxation rate Pp. 

The momentum scattering time Tp(n) is expected to 
decrease with increasing carrier density n. Therefore, the 
spin relaxation time, r^, should increase with increasing 
n if spin relaxation is dominated by the DP mechanism 
and decrease with increasing n if spin relaxation is dom¬ 
inated by the EY mechanism. In our experiments, at 
T= 10 K, we found the initial decay of Pc is character¬ 
ized by Tg oc n“® ®®. The decreasing spin relaxation time 
with the increasing density (i.e., decreasing Tp) suggests 
that the EY mechanism plays a larger role than the DP 
mechanism in the spin relaxation of cold excitons near 
the band edge in GaSe at low temperature. 
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V. THEORETICAL MODELING 


The dynamics of resonantly excited nonthermal ex- 
citons in quasi-two-dimensional systems such as GaAs- 
based quantum-well structures are affected by several 
physical processes^®’'^®^®®: (1) momentum relaxation of 
excitons, (2) spin relaxation of excitons, and (3) the en¬ 
hanced radiative recombination and propagation of ex- 
citon polaritons. For nonresonantly photoexcited carri¬ 
ers, one should also consider the contribution to dynam¬ 
ics from free carriers, particularly at high temperature. 
Here, we calculate the population of excitons in various 
spin and momentum states in GaSe as described in Sec. 
II using a simplified exciton spin-flip model (Fig. 5). The 
model is adapted from a unified model for resonantly ex¬ 
cited excitons in GaAs-based quantum wells®^. In prin¬ 
ciple, this model should be valid only for resonant pho¬ 
toexcitation at low temperature where the contributions 
from free carriers is negligible. Nevertheless, we find that 
this model reproduces the essential experimental photo¬ 
luminescence polarization properties and dynamics. 

We label exciton states as |±) = |1,±1), |1,0), and 
|0,0). We first consider the case when the excitons are 
photoexcited in non-radiative, high-momentum states 
(ATy > Kq) labelled, for example, \+k), where ATy is 
the in-plane momentum of excitons and Kq is the pho¬ 
ton momentum. Exciton spin-flip (with rate Wx) trans¬ 
fers population between |-|-) and |—) states, while elec¬ 


tron/hole spin-flip (indistinguishable in our experiments 
and so characterized by a single rate (Wg)) populates the 
dipole-inactive |1,0) non-radiative (dark) state, and the 
singlet |0,0) state (dipole-active for E || c). Following 
Vinattieri et al.®^, we divide the manifold of ATy states 
into two sets, one for nearly zero ATy and the other for 
finite large ATy states. Each set includes the four exci¬ 
ton states. Absorption and emission of acoustic phonons 
induce transitions between these two sets. We consider 
only spin-conserving transitions with an effective scat¬ 
tering rate W^- To simplify the model, we neglect any 
thermal factors associated with the spin-flip rates of ex¬ 
citons and electron/hole. 

The time-dependent population in each state is given 
by a set of coupled equations: 


—W — Mij Nj -|- G{t) S+ik,i 


( 5 ) 


where W is the column vector 
(Ak|_i, 7V_i, Wo, Aoo, A_|_ifc, 7V_ifc, Wofcj Aoofc) and M is 
a 8 X 8 matrix. N^i, N_i, Niq, Nqq are populations of 
K\\ < ATo states |-|-), |—), |1,0), and |0, 0), respectively. 
iV+ife, A_ifc,7Vio/c, A^oofc are corresponding ATy > ATq 
states. M is the following matrix: 


M = 


A C 
D B 


( 6 ) 


where A, B, C, and D are the following 4x4 matrices: 


A = 


-{WR+Wx+Ws+Wkp 

Wx 

Ws/2 

WJ2 


Wx Ws/A 

-{WR + Wx+Ws + Wkp) W,/i 

Ws/2. -{Wj2+Wkp 

Ws/2 0 


Ws/A 

WsjA 

0 

-{Wl,+W,l2+Wkp) 


B = 


- -(Wpp+Wx+W,+Wkm.) 
Wx 
WJ2 
Ws/2 


Wx Ws/A Ws/A 

-{W„r- + Wx+Ws + Wkm) Ws/A Ws/A 

Ws/2. -{Wpr- + Ws/2 + Wkrp) 0 

Ws/2 0 -(Wpp+Ws/2+Wkrp)- 


^^km 

0 

0 

0 ■ 


[Wkp 

0 

0 

0 ■ 

0 

Wkm 

0 

0 

,D = 

0 

Wkp 

0 

0 

0 

0 


0 

0 

0 

Wkp 

0 

0 

0 

0 

Wkm. 


0 

0 

0 

Wkp. 


Wfi and are the radiative recombination rates for |±) 
and |0,0), respectively. Wx and Wg are the exciton and 
electron/hole spin-relaxation rates. The acoustic-phonon 
scattering rates Wkp and Wkm are defined by Vinattieri 
et al.®® as 


Wkp = Wk exp 


hVh 

keT 


Wkm = Wk 


^1 — exp 


keT 


( 7 ) 


where Wk is the effective scattering rate with phonons, 
and hTh is the homogeneous linewidth. The spin-flip 
rates of electrons and holes cannot be distinguished in 
this model because spin-flip of electron and holes results 
in identical transitions within the model. 

There are six parameters in the model (VF)j, Wn, Wx, 
Ws, Wk, Th). However, the present and prior experimen¬ 
tal data strongly constrain the values these rates can take 
in fitting the experimental polarized PL dynamics with 
this model (Fig. 5). The radiative recombination rate 
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FIG. 5. Schematic of the model for the exciton dynamics. 
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FIG. 6. Calculated polarized PL dynamics for T = 10 K and 
djL = 650 nm. 


for state |0,0) is set to 30 Wr based on the relative 
absorption coefficients between E || c and E _L c light^^. 
The value of is not available but in principle can 
be measured independently from temperature-dependent 
PL linewidths®'^^®^. For simplicity, we set HTh = ksT 
because the linewidth increases from ~30 meV at T = 
10 K to ^50 meV at T = 300 K. Wr is then largely de¬ 
termined by the decay of population (the measured PL 
decay). 

We first consider the total time-dependent PL un¬ 
der linearly and circularly polarized excitation condi¬ 
tions. Under linearly polarized excitation, the photoex- 
cited electrons and holes are assumed to be equally dis¬ 
tributed over their respective spin states, i.e. |±l/ 2 ), 
during the initial 2-ps laser excitation. The nongemi¬ 
nate (bimolecular) formation of excitons then produces 
an initial exciton population distributed equally over 
the three triplet exciton states |±) and |1,0). We con¬ 
trast this to the case of circularly polarized excitation 


(tT+), under which only the bright exciton state |-|-) 
is initially formed. The total time-dependent popu¬ 
lations of |-|-) and |—) created by circularly and lin¬ 
early polarized light are labeled I^{t) and re¬ 

spectively. Then the ratio R{t) = /+(t)/J^(t) = 
[/■'■/■'■ (<) -I- /+/“(<)] / -I- will decrease 

from approximately 1.5 to 1 as the populations in 
the three triplet states eventually become nearly equal 
through spin relaxation under cr'*' excitation. 


In order to estimate Wx-, we first focus on polarized 
PL dynamics under circularly polarized excitation. The 
sum of time-dependent cr"'' and a~ PL under circularly 
polarized excitation and cr^ and PL under linearly 
polarized excitation are both independent of Wx 

because Wx does not change the total population in the 
two radiative states. By fitting experimental results for 
the polarized and total (i.e., unpolarized) PL versus time, 
we find that polarized PL dynamics is dominated by Ws 
at cryogenic temperatures (T = 10 K). 


The preceding discussion of the decay of spin polar¬ 
ization applies to band-edge carriers. However, it does 
not account for spin relaxation of carriers excited above 
the gap, which exhibit faster spin relaxation as they cool 
and scatter into the region iLy < Kq^^. At T = 10 K, 
the time-dependent degree of circular polarization {pc{t)) 
exhibits biexponential decay. The corresponding time- 
dependent PL under linearly polarized excitation also de¬ 
cays biexponentially. In addition, the ratio R(t) decays 
to 1 within about 20 ps. We can account for these exper¬ 
imental features without appeal to the detailed nonequi¬ 
librium carrier distributions, by assigning a spin-flip rate 
that decreases linearly from about W^ « 500 Ws to Ws 
during the thermalization period Tth- The thermaliza- 
tion time can be determined by fitting to the polarized 
PL dynamics at low temperature [rth ~ 80 ps at T = 10 
K). After including such phenomenological parameters 
Wg = 500 Ws and Tth = 80 ps in the model, we reproduce 
quantitatively Pc{t): and and 

qualitatively [R{t) — 1] (Fig. 6). The parameters used 
for the calculations are as follows: Wr = 0.006 ps“^, 
Wnr < lO-'^ pS-\ Wx < 10-^ pS-\ Ws = 0.001 pS"!, 
Wg = 0.5 ps“^, Tth = 80 ps, HTh = 2 meV, and Wk = 1.0 
ps“^. The calculated PL dynamics (Fig. 6) agree with 
the experimental results (Fig. 4). 


This model is also consistent with the remote-edge- 
luminescence (REL) dynamics reported previously^®. 
Excitation with E T c light results in linearly polar¬ 
ized emission with polarization along the line from the 
excitation spot to the emission spot. The highly linearly 
polarized REL is due to scattering of excitons from the 
|±) excitons into the |0, 0) state. The rise of the REL 
corresponds to the optical spin decay of the PL from the 
focal spot (with a time constant ~20-30 ps). 
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VI. SUMMARY 

We have analyzed exciton spin dynamics in GaSe under 
nonresonant circularly polarized optical pumping with a 
phenomenological exciton spin-flip rate-equation model. 
The model reproduces the initial unity circular polariza¬ 
tion as well as the biexponential decay with a sub-20 ps 
and a >500 ps time constant under nonresonant optical 
pumping with excess energy up to 0.15 eV at T = 10 K. 
The separation of the non-degenerate conduction and va¬ 
lence bands from other bands results in angular momen¬ 
tum preservation for both electrons and holes. Spin re¬ 
laxation of hot spin-polarized carriers nonresonantly op¬ 
tically injected into noncentrosymmetric e-GaSe is dom¬ 
inated by the DP spin-relaxation mechanism; a sub-20 


ps spin-relaxation time results. By contrast, the spin- 
relaxation rate of cold excitons {e-h pairs) formed near 
the T-point near the band edge is greatly reduced as a re¬ 
sult of the suppressed DP and EY spin-relaxation mech¬ 
anisms. 
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